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Abstract. The theory of "subalgebra basis" analogous to standard basis (the 
generaUzation of Grobner bases to monomial ordering which are not necessarily 
well ordering [l].) for ideals in polynomial rings over a field is developed. We 
call these bases "SASBI Basis" for "Subalgebra Analogue to Standard Basis 
for Ideals" . The case of global orderings, here they are called "SAGBI Basis" 
for "Subalgebra Analogue to Grobner Basis for Ideals", is treated in Sasbi 
bases may be infinite. In this paper we consider subalgebras admitting a finite 
Sasbi basis and give algorithms to compute them. The algorithms have been 
implemented as a library for the computer algebra system SINGULAR [2]. 



1. Introduction and Preliminaries 

Let K he a, field and K[xi, . . . , Xn] the polynomial the ring over the field K in n 
variables and . . . ,Xn]] the formal power series ring. Let G C A^\{0} where 

Ai is the maximal ideal of iir[[a;i, . . . , a;„]]. We define 

K\[G]] = {(3(51, ...,g,)\Qe K[[yi, . . .,ys]] and gi,...,gs eG for somes}. 

In analogy to the theory of Grobner bases for ideals in K[xi,. . . ,Xn] resp. stan- 
dard bases for ideals in /4r[[xi, . . . ,Xn]] there exist a theory of Grobner bases for 
subalgebras of type K[G] called Sagbi basis (c.f [6]) and of standard bases of sub- 
algebras of type if[[G]] developed in [3]. Let fi, . . . ,fm £ K[xi, . . . , Xn] and con- 
sider the ideal / — {fi, . . . , fm)K[[xi, . . . ,Xn]]- Suppose we want to compute a 
standard basis of /. There are at least three possibilities. Using Buchberger's al- 
gorithm for well orderings we can compute it upto a given degree. There exist a 
theory of standard bases in K[xi, . . . , Xn](xi,...,x„) induced by a local degree or- 
dering (see [1]) and we can compute a standard basis gi, . . . ,gs € K[xi, . . . ,Xn] 
of (/i, . . . , fm)K[xi, . . . , Xn]i^xi,...,x„) using Mora's tangent cone algorithm (c.f [1], 
[5]). It can also be computed via homogenization (Lazard's algorithm c.f j^). In 
this paper we will developed a subalgebra analogue of standard bases of ideals in 
K[xi,. . . , Xn](xi,...,x„) ■ We will introduce subalgebra bases in a suitable localization 
of [G], which we call it Sasbi bases. They can be computed upto a certain degree 
similar to the standard basis case. They can be computed via a Sagbi bases if the 
homogenized algebra has a finite Sagbi basis. In this case they can also be com- 
puted directly using a generalization of the tangent cone algorithm. We will show 
that Sasbi bases are subalgebra standard bases of fir[[G]]. Using homogenization 
we will get a finiteness condition. If i^[G''] has a finite Sagbi basis then i^[[G]] has 
a finite subalgebra standard basis. The aim of this paper is to give an algorithm to 
compute these bases provided G C i4'[[a;i, . . . , Xn\\ is a finite set and the subalgebra 
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standard basis for i^[[G]] is finitcQ- We use the notations from [1] and repeat them 
for the convenience of reader. 

Definition 1.1. A monomial ordering is a total ordering > on the set of monomials 
Moun — {x" I a £ N"} in n variables satisfying 

for all a, /3, 7 6 N". We also say > is a monomial ordering on ii'[a;i, . . . , a;„] meaning 
that > is a monomial ordering on Mort„. 

Definition 1.2. Let > be a fixed monomial ordering. Write / G if [xi, . . . , 
/ ^ 0, in a unique way as a sum of non-zero terms 

/ = acx" + fl/jx^ + . . . + a^x'^ , x" > a;'^ > . . . > x^, 

and Oq, . . . a-y G if. We define: 

1. LM( f) := x", the leading monomial of /, 

2. LE{f) := a, the leading exponent of /, 

3. LT{f) :— Gax", the leading term of /, 

4. LC{f) :— tta, the leading monomial of /, 

5. tail{f) := / - LT{f). 

6. ecartif) := deg{f) - deg{LM{f)). 

7. support{f) {x" , x^ , . . . , x^}, the set of all monomials of / with non-zero 
coefficent. 

8. ordif) - deg{LM{f)). 

Definition 1.3. Let > be a monomial ordering on MoUn- 

1. > is called global ordering if a;" > 1 for all a 7^ (0, . . . , 0). 

2. > is called local ordering if < 1 for all a ^ (0, ... , 0). 

3. > is called local degree ordering, if > is a local ordering and 

a;" > x'^ ^ deg{x°') < deg(x'^) 

Definition 1.4. Let M G GL{n, R). We can use M to obtain a monomial ordering 
by setting 

x" >A/ x'^ ■.<=^ Ma > MP, 
where > on the right-hand side is the lexicographical ordering on R". 

Theorem 1.5. (c.f. [1], page 18 )Any monomial ordering can be defined as >m 
by a matrix M £ GL(n,R). 

Definition 1.6. For any monomial ordering > on Mon{xi, . . . ,x„), we define a 
multiplicatively closed set. 

5> := {u e K[xi, . . . , x„]\{0} I LM{u) ^ 1} 

Let K[xi, . . . ,x„]> := S>^^K[xi, . . . , x„] = | /, u G K[xi, . . . ,x„],LM(m) = 1} 
the localization of K[xi, . . . , x„] with respect to S'> and call iir[xi, . . . , x„]> the 
ring associated to iir[xi, . . . , x„] and >. 



They may be infinite as /^[[a;!, . . . , x^]] doesn't satisfy the ascending chain condition with 
respect to subalgebras. 
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Definition 1.7. Let > be any monomial ordering. For / G K[xi, . . . ,Xn]> choose 
u G K[xi, . . . , Xn] such that LM{u) = 1 and uf G K[xi^ . . . , Xn]- We define 

LMif) = LM{uf) 

LC{f) = LC{uf) 
LT{f) = LT{uf) 
LE{f) = LE{uf) 
tail{f) = f- LT{f) 

Definition 1.8. Let G be a subset of K[xi, . . . , Xn] 

A G-monomial is a finite power product of the form G" = g"^ . . . g^™ where gi € G 
for i = 1, . . . , m, and a = (ai, . . . , am) S N™. The set of all G-monomial is denoted 
by: 

Mona = {G"| a e N™, m G N} 
Let G = {gi, . . . , gs} C K[xi, . . . , Xn] and > be a local ordering. We define 
K[G]> = (5> n if [<?!, . . .,gs])-'K[g,, . . . ,5,]. 

2. Sasbi basis of if [[G]] 

Let tVJ is a maximal ideal of if [[xi, . . . ,Xn]]- We fix a local degree ordering > 
and use the notation of definition 1.2 which make sense in if [[xi, . . . ,Xn]] too. In 
this section we recall some result of [3] and give an algorithm which computes Sasbi 
basis if [[G]] upto a certain certain degree. 

Definition 2.1. Given two elements g,h G K[[xi, . . . ,Xn]], we will say that g 
reduces to h with respect to G if their exist G-monomial G" and 7 G if such that 

h^g- 7G", with h = OoT LM{h) < LM{g) 

In this case we will write 

g^h, 

and we have that g — /) G if [[G]]. 

Consider a chain (possibly infinite) of reductions 

g ^ ill ^ /l2 ^ ■ ■ • /im ^ • • ■ 

This implies there exist G-products G"' ' and a^a) G if \{0} such that 

m 

hm ^ g — ^ flaCi) G" , 
1=1 

and because of the definition of the reduction 

LM(a„(i)G""') > LM(a„(2)G"*'') > . . . 
If the chain is infinite, we get the following sequence in if [[X]] : 

m 

Sm = ^a„(i)G" ' , m > 1. 
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This sequence happens to be convergent in iir[[xi, . . . with respect to the 
A^-adic topology. We denote the Hmit of the sequence {sm)m>i by s. Since aU the 
terms are in the complete subalgebra ^^[[0]] so we have that s G iir[[G]]. 

Definition 2.2. If the reduction ft, of g is zero or for all e support{h),x^ ^ 
LM{G") for all G-monomial G" then h is called Normal form of g with respect 
to G. We denote the normal form h by NF{g \ G). 

Normal forms always exist but the computation may use infinite reductions. For 
computational reason we give an algorithm which computes the normal form up to 
the degree d. 

Algorithm 2.3. Let > be any local degree ordering in K[[xi, . . . 
Input: G C M\{0} , g £ K[[xi, . . ■,Xn]], deZ. 

Output: h—NF(g \ G, d) (the normal form of g with respect to G up to degree dj^. 

• /i := g; 

• while{h 7^ and ordih) < d) 

n - {G"| LAf(G") = LM{h)} ^ <f>; 
ifTh^(i> 

choose G" e Th ; 

h — h — LC{h) . 

else 

return {LT{h) + NF{h - LT{h) \ G, d); 

• end(while) 

• return h; 

Definition 2.4. We say a set G C M\{0} is Sasbibasi^ of K[[G]] if 

K[L{K[G])]^K[LM{G)] 

where L{K[G]) = {LM{g)\g e K[[G]]\{0}}, i.e G that is a sasbi basis if for aU 
/ G K[[G]]\{0}, 

LM{f) = LM{G°') 

for some G-monomial G". 

Example 2.5. The set G = {a^^, X^i^s C i4r[[a:i, . . . , is a sasbi basis for 
K[[G]]. Indeed, if g e K[[G]]\{0}, then LM{g) = 1 or LM{g) = x", for some 
a > 2. Hence LM{g) G K[x^,x'^] = K[LM{G)]. 

Now for the characterization of the sasbi bases in iir[[xi, . . . , x„]] similar to those 
in K[xi, . . . , Xn] we need to define an analogue of the S-polynoniial. 

Definition 2.6. Let G — {gi, . . . ,gs} C A1\{0}. An S-polynomial is an element 
of the form 

aG" ~bG^ 

where a,b E K\{0} and G" and G^ are G-monomials, such that LT{aG°') — 
LM{bGl^). 

^for theoretical reasons we allow G to be infinite and d = oo. We have seen that for d ^ oo 
the normal form NF{g \ G, d) converges in the A^-adic topology. We call this limit NF{g \ G), the 
normal form of / with respect to G. 

''in [3] this is called a standard basis of subalgebras. We use this notation to be similar to 
sagbi bases introduced in [B]. 
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Next theorem gives criteria for a set to be a sasbi basis of iir[[G]]. 

Theorem 2.7. (c.f. [3], page 50) Given G — {51, . . . ,gs} C A^\{0}, G is a sasbi 
basis of K[[G]] if and only if every S-polynomial of G has a vanishing normal form 
with respect to G. 

The following is the analogue of Buchberger's Algorithm for subalgebras in 

K[[xi, . . .,Xn]]- 

Algorithm 2.8. 

Input: A finite subset G of Ai\{0}. 
Output: A sasbi bases F for K[[G\]. 

• F = G ; 

• oldF = (j) ; 

• while (F ^ oldF ) 

S := {s\s is a S-polynomial ofF}; 
R:={r\r = NF{s € S \ F) andr ^ 0}; 
oldF = F; 
F = FUR; 

• return F; 

3. Sasbi Basis in the Localization of K[xi, . . . ,Xn] 

In this section first of all we will introduce Sasbi bases in _fir[G]> and prove that 
Sasbi bases in iir[G]> are also Sasbi bases in if[[G']]. We will prove that also in the 
general case the computation of a Sasbi basis with respect to a local ordering can be 
reduced using homogenization to the computation of a Sagbi basis with respect to a 
suitable global ordering. This is also here a very expensive way to compute a Sasbi 
basis. Therefore later a more efficient algorithm is presented. We introduce notion 
of weak sasbi normal form of a polynomial with respect to G in K[xi, . . . , a;„]> and 
give an algorithm to compute it. Then we give a criterion for a set to be a Sasbi 
basis, which is the base of an algorithm to compute the Sasbi basis. 

Let G = C K[xi,...,Xnl G = K[G]y and L{g) = {LM{g)\g e 

G\m}- 

Definition 3.1. A subset S C G called SASB^^ Basis of K[G]> if 

K[LiG)] = K[LiS)] 

i.e for aU g e i^[G]>\{0} 

LM{g) = LM{S°') 

for some S-monomial S*". 

If > is global, a Sasbi basis is also called a Sagbi basis. 

If we just say that S* is a Sasbi basis, we mean that S" is a Sasbi basis of the /4r[S']> 
generated by S. 

Theorem 3.2. Let K[xi, . . . , Xn]> C . . . , Xn]] be equipped with local degree 

ordering >. Let G = {.gi, . . . ,gs} a subset of K[xi, . . . ,a;„]. If S is a Sasbi basis 
of K[G]y then S is a Sasbi basis of K[[G]]. 



SASBI stands for "Subalgebra Analogue to Standard Basis For Ideal" 
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Proof. For g e we have to prove that there exist G-nionomial G" such 

that LM{g)=LM{G"). If g £ K[[G]] there exist H e K[[yi, . . . .y^]] such that 
we have g = H{gi, . . . ,gs)- There exist a decomposition oi H = H^^^ + H^^\ 
ir(°) e K[yi, ...,ys] and i/fi) G K[[yi, . . . such that 

LMiH^°\g,, . . .,g,)) = LM{H{g„ . . .,g,)) = LM{g). 

Since S is Sasbi bases of K[G\y there exist a G-monomial G", such that LM{G°') = 
LM{H^°\gi, . . .,gs)). We get LM{g) = LM{G") which shows that 5 is a Sasbi 
bases for ^^[[G]]. □ 

Now we want to show how to reduce the SASBI bases computation for local 
orderings using homogenization with respect to a variable "t" to the computation 
of SASBI bases for global orderings. 

Theorem 3.3. Let H ^ {gi, . . .,gm} C K[xi, . . . ,a;„] and K[H] = K[gi, . . .,(/„]>. 
Here > is a local monomial ordering given by a matrix M. Consider K[t, x\,. . . , a;„] 

/II... 1 \ 


with monomial ordering >h defined by the matrix 

: M 
\ / 

>h is a global ordering. We define Gi to be Gi := g^ € K[t,xi, . . . ,Xn]. Assume 
S = {Si,S2, . . . , Sk} C K[t, xi,X2, ■ ■ ■ , Xn] is a Sagbi basis of K[Gi, . . . , Gm] with 
respect to >h. Let Sj := Sj{t = 1) , 1 < i < fc, then S = {si, . . . ,Sfe} is a Sasbi 
basis of K[gi,g2, . ■ ■,gm]>- 

Proof. We want to show 5 is a Sasbi basis for K[H]y. For this we have to show 
that 

1. 5* C K[H]y. 

2. For g e K[H]y there exist a = {ai,a2, ■ ■ ■ ,ak) G N*" such that LM{g) = 
LM{S°'). 

1) We know that S = {Si, . . . , Sk} is a Sagbi basis of K[Gi, . . . ,G„i] so Si = 

jG"''J with 7j j s K. Put t = 1 we get Sj = j^"''' this implies Sj € 
K[H]>- 

2) For g G -ftr[if]> there exists u e S> D K[H] such that u.g = then 
there exists p G Z such that tP.u^.g'^ = Y.liiH'^'Y = T.ljG"^ . We have that 
S is Sagbi basis of K[Gi, . . . , Gm]- Then there exists a — {ai,a2, ■ . ■ , Qk) G N'' 
such that LM{tP.u''.g'') = LM{S"). Since LM(G)|t=i = LM(G|t=i), therefore 
LM{tP.u''.g'')\t=i = LM{g), since LM{u) = 1 as u e S> and LM{S")t=i = 
LM{S"), we obtain LM{g) = LM{S"). □ 

Theorem 3.3 shows that Sasbi bases are computable in many cases. It turns out 
that similar to the theory of standard bases with respect to local orderings for ideals 

this approach is not very efficient. Therefore one should like to have an efficient 
way for computing Sasbi bases. The basis for this is the concept of the normal 
form. 

Definition 3.4. Let G and g be a finite subset and a polynomial in K[xi, . . . , x„] 
respectively,such that _fi'[G]> admits a finite Sasbi bases and g G /^[G]. 
We say that a polynomial /i is a Weak SASBI normal form of g with respect to 
G, and we write h = SNF{g\G), if 
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0. h = SNF{0\ G) = 

1. /i ^ ^ LM{h) i K[LM{G)] 

2. There exist unit u € riK[G] such that ug — h ha-s a representation with respect 
to G, that is either ug — h = or ug — h = J2i=i 7*^"' where ji € K and LM{g) = 
max-Lj{LM(7iG"*)}. This representation is called SASBI representation. 

Algorithm 3.5. 

Input: f , G, > a local m,onomial ordering. We assume that G — {gi, . . . ,gs} and 
f are subset and polynomial in K[xi, . . . , x„] such that f G K[G]. We also assume 
there exist a finite Sagbi basis of K[H] where H = G^ the homogenization of G 

with respect to "t" , a new variable. 

Output: h ^ K[x] a polynomial weak Sasbi normal form of f with respect to G. 

• h:=f 

• T:^G 

• while{h ^ and Th = {T", T-monomial\ LM{T°') = LM{h)} ^ </> 

choose T" e T/j such that ecart{T°') is minimal; 
ifecart{T") > ecart{h) 

T:=TU{/i}; 
h = h- LC(Hl) T°'for some^i e K; 

• return h; 

Proof. Termination is most easily seen by using homogenization: start with h := f^ 
andH := G^ = {g^\g gG}. 
The while loop looks as follows 

• while(/i ^ and Th = {H", iJ-monomial | LM{H") = tl^LM{h) for some /?}) 

choose g e T/, in a way with /? > is minimal; 
if /3 > 

T = TU{h}; 

h := (/7.k=i)^ 

By our assumption K[H] has a finite sagbi bases, there exists some positive integer 
A'' such that K[L{Ty)\ becomes stable for v > N, where Ty denotes the set T 
after the v-th turn of the while loop. The next h, satisfies LM{h) E K[L{Tn)] = 
K[L{H)], whence LM{h) = LM{H°') for some if" G K[H] and /? = 0, that is, r„ 
itself becomes stable for v > N and the algorithm continues with fixed T. Then it 
terminates, since > is a well ordering on K[t, x\. 

To see the correctness, consider the i-th turn in the while loop of algorithm. 
There we create Tj = {ffi,g2, ■■■,ga, ho, h\, . . . , hi-2} such that hi = hi-i — 7*T" 
and LM(T«''-') = LM{hi-i) > LM{h,) where T"*'' is T,- monomial. 
Suppose, by induction, that in the first i — 1 steps we have constructed SASBI 
representations 

Ujf = J2 li'^G"''^ + hj where 7/ G K and LM{f) = maxJ'fi{LM(7PG«''')}. 

1=1 

where uj e S> n K[G] and 1 < j < i - 1 

We have to prove 3ui e S> n K[G] and Uif = 7^ ^G"" + hi and LM{f) = 

1=1 
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maxfi:i{LM(7f^G"J")}. 

We have two possibilities 

1) T"^" = G"'"' is G-monomial. 

2) T«"' = 17''' is Tj-monomial. 

Induction step: Consider the SASBI representation for j = i — 1. 

Ui-if= ^ 7f'''G"'''"''+/i«-i and LM(/) = max,tT" {iM(7f -^^G"'*"'" )}. 
For the first case in induction step , replace hi-i by 7'G"''^ + hi ,and obtain 

= "^'7(-i)G«r%yG«'%/.,. 

Put = and 7f-^'G"^'' = 7«G"J*' 1<1< v^'-^\ 7^G"''' = 7i:^G"^'^' we 
get the required representation 

1=1 

As LM(7^fG"") < LM(7f^G""), 1 < / < - 1, from this condition we get 

LM{f) = niax|!i'J{LM(7}'^G"! )} which shows representation is Sasbi. 

For the second case in induction step replace the hi-\ by 7'T'" + /i^ , it becomes 

1=1 

We can write Tf^ = G^^ where H = {/lo, /ii, • • • , /ii-2}. Since we are in 

second case not all the components of 7^*-' are zero. Since LM{hi^i) < LM{hj) 
ioY j <i-2 and LM{T'^'-'^) = LM{h,_i) < LM{f) it follows that LM{GP^^ < 1. 

Since hj = Ujf — J2 Ji'^^G"'^ ,1 < j < » — 2 , we can replace h',s by this expression 
1=1 

therefore 

Ti''^'^ = G^Wi?(uo, Ml, ... , Ui_2, /, ffi, . . . , 9s) f + G^^^L-^^'K 

For a suitable polynomial and L = { E if^G"'" , • • • , E ll'"^^^"'' ' }• 

Since uq, ui, . . . , Ui-2, / e i^fG] and LM{Gl^(^) < 1 it follows that 

= - fG'^'^^R e 5> n iffG]. 

Since LM(7;^''^G"'^' ) < LM{f) it follows that leading monomial of any G-monomial 
occuring in G'^^'^L'''*'*) is smaller than the leading monomial of /. This implies 

„(-i) 

uj= J2 7r'^G«^%yG«'''L^W+/ii. 
1=1 

is a sasbi representation since LM{f) = max-^i'i ^'{LM(7|'~^^G"' ')}. □ 
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Example 3.6. In the localization of the univariate polynomial ring K[x\y where 
> is the local ordering take g — + a;^ and G ~ {gi — x^ + x^ , g2 = x — x^} we 
want to compute the weak Sasbi normal form of g with respect to G. 
In the first reduction we select the G-monomial gi = x^ + x^ with minimal ecart 
such that LM{ho) — LM{f) =^ LM{g) ^ x^ , we have ecart(hQ)=l, ecart{gi) = 3, 
so ecart{gi) >ecart(ho) therefore we have to enlarge G — {gi = x'^ + x^,g2 — 
X — x^ , g3 = x^ + cc'*} and 

hi =hQ~ gi 

In the second reduction we select the G-product = {x — x'^){x^ + x'^) = 

x^ — x^ with minimal ecart such that LM{hi) = LM{g2g3) = x'^ . Now we have 
ecart{g2g3) = 2, ecart{hi) — 2 so G remains the same and 

h2 = hi - 92.93, 
we get h2 — 0. Now we summarize and obtain 

/l2 = /ii - .9253, 

= x^ -x"^ -(x-x^){x^ +x^). 

As hi ^ g ~ gi 

h2= 9~ 9i~ 9293, 
= + - (x^ - x^) -{x- x^){x'^ + x^), 

93^9 we get 

9-992 = 9i + ^2, 
x^ + x^-{x- x^){x^ + x'^) ^x^ - x^ 

(1 - 52)5 =91+^2 

{l-x-x^){x^+x^) = x^~x^. 

we have 1 — g2 = 1 — x — x^ Sy Cl K[G] and h2 = is the weak Sasbi normal 
form. 

Definition 3.7. Let G ~ {gi, . . . ,gm} C K[xi, . . . ,Xn]- Let 

ARiG) := {h e K[yi, y™] | h{LM{gi), LMig„,)) = 0} C K[yi, y,„]. 

Definition 3.8. Let G C K[xi,...,Xn] and EjLi7»G"' G K[G]. We define the 
height ht{J2'^^ll^G°'') = max'^^j^{LM{G°'')}. 

Theorem 3.9. (SASBI basis criterion) Let G={gi, g2, ■ ■ ■ ,.9m} be a subset of 
K[xi, ...,Xn]. Assume that K[II] has a finite sagbi basis where H = {gi^ , ■ ■ ■ ,gm^}, 
the homogenization ofG. LetS {-Pi, Pk) be a generating se0 of AR{G) . Then 
G is a SASBI basis for K[G]> if and only if for each 1 < j < k, SNF(Pj (G)\ G)=0. 



'The set of 5-polynomials defined in definition 2.6 defines a generating set of AR{G). 
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Proof. (^) Suppose that S'A^i^(Pj(G) \ G) ^ 0. This implies that LM(5iVF(Pj (G)| 
G) ^ K[LM{G)] by the property of the weak Sasbi normal form . We have 
Pj{G) G K[G] therefore SNF{Pj{G) \ G) £ K[G]. Since G is a SASBI basis 
of K[G]> we have LM{SNF{Pj{G) \ G) e K[LM{G)]. This is contradiction to 
the assumption that LM{SNF{Pj{G) \ G) ^ K\LM{G)\. 

(<^=) To prove that G is SASBI basis , we have to prove that g S K\G\-^ has a 
SASBI representation with respect to G, that is there exist u G S'> fl iir[G] such 
that 

ug = ^7iG"* with LM{g) = ht{J2l^G"') 

i=l i=l 

Let g G if [G]>, choose u G Sy r\K[G] such that ug = Y^i^i jiG""' , furthermore, 
we assume that this representation has the smallest possible height of all possible 
representations of ug in K[G]. We denote this height by X:=max'-^i{LM{G°'')}. 
It is clear that LM{g) < X. Suppose that LM{g) ^ X. Without loss of gener- 
ality, let the first jj, summands in the above representation of g, be the ones for 
which X=LM(G"'). Then cancelation of their leading terms must occur, that is, 
'^i-i%LT{G°'') = 0, and hence we obtain a polynomial in K[yi, ...,y„i], P{y) = 
liU"' G AR{G). Since, S — {Pi, ...,-Pfe} is a generating set of AR{G) we can 

write 

k 

For suitable fj G K[yi, . . . , ym]- Furthermore , note that 

ht{P{G)) = max']^iht{gj{G))ht{Pj{G)) = X 

where, fj{G) and Pj{G) are considered as expressions in the gi 's 
On the other hand : 

By assumption we have for all 1 < i < fc , SNF{Pj{G) | G) = 0, which means that 
WjPj{G) has a SASBI representation, WjPj{G) = Y^'iLiUjG"'^ , for suitable Wj G 
S>nK[G] and LM{Pj{G)) = max^'^ JLM(G"'0} S ht{Pj{G)). The inequality is 
strict since Pj G AR{G), we may assume that w = wj, where 1 < j < fc . For each 
j, we have 

wfj{G)PjiG) = Y^liMG)Gl. {**) 

i=l 

If we define Xj to be the height of the right hand side in the equation , then 
obtain 

Xj $ max']^iht{fj{G).ht{Pj{G)) = X. 
Finally, the equations (*) and (**) imply that : 



V 

ug = P{G)+ ^ 7iG' 
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i=i 1=1 i=n+i 

^ V ' ^ V ' 

sumi sum2 

If we examine the expressions of the above equation, we see that Xj < X ; for all 1 < 
j < k therefore ht(sumi) = max^^^Xj < X. By the choice of /i, ht{sum2 < X). 
But this contradicts our assumption that we have chosen a representation of h with 
smallest possible height. Thus, G is a SASBI basis of K[G]>. □ 

This theorem is the base of following algorithm : 

Algorithm 3.10. Let > be a local monomial ordering on K[xi, 

Input: A finite subset G C K[xi, Xn]. Assume K[G]y admits a finite SASBI 
basis and K[H] admits a finite sagbi basis where H = G^ is the homogenization 
with respect to new variable "t". 
Output: A SASBI basis F for K[G]>. 



• F = G; 

• oldF = (j); 

• While (Fy^oldF) 

Compute a generating set S for AR{F); 
V = S{F); 

Red={SNF(p \F)\p&V\{Q}}; 
oldF = F; 
F = FU Red; 



• return F; 

Example 3.11. LetG = {gi = x/,g2 — x4:+x5+x6,g3 — = a:'', 35 = y^+x^} 

is a subset K[x,y] and > the degree lexicographical local monomial ordering. We 
consider K[G]y = K[x^,xA + xb + x^ , y"^ , x"^ , + a;^]> . Then we have an ideal 
AR{G){G) — (si — x^y^ + si — x^ -\-x^). We can take the reduction of si by 
91^95 (with minimal ecart) we obtain 

h = si- gi^gr, 

xV + \x''-{x'f{y^ + x') = Q 

so SNF{si I G) = 0. There is no G-monomial G" such that LM{G") = LM{s2) = 
x^ so SNF{s2 \G) = x^ + x^. We have new G = G U {ge ^ x^ + x^}. Then we 
have new AR{G){G) = (si = x^y'^ + ^x^^,si = x^ + x^), so obviously weak sasbi 
normal form of si and S2 are 0. This shows that G = {91, g2, 93, 94, 95,96} is a 
Sasbi basis of K[x^, xA + + x^,y'^, x^ , y^ + a;^]> . 

We have presented the theory of sasbi basis for /ir[G']>, where G is finite subset 
Xi, . . . ,Xn] and > is local orderings, but it is still an open problem for mixed 
orderings. 

4. Implementation in SINGULAR 

In this section we will give an overview of the main procedures which we have 
implemented in SINGULAR. In this overview we will present these procedures and 
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give by concrete SINGULAR examples to explain their usage. We have imple- 
mented three types of procedures: 

1) Weak sasbi Normal form procedure 



" WSNF procedure" : It is an implementation of Algorithm 3 (ecart driven nor- 
mal form) to obtain weak sasbi normal form of a polynomial. 

SINGULAR Procedure: 

LIB"algebra. lib" ;// we need this library for "algebra.containment" 
// procedure 

proc WSNF (poly f .ideal I) 

{ 

ideal G=I ; 
poly h=f ; 
poly hl,j ; 
list L ; 
map psi ; 

while (h!=0 && hl!=h) 
{ 

L= algebra_containment(lead(h) ,lead(G) ,1) ; 
if (L[l]==l) 

{ 

def s= L[2] ; 

psi= s,meLxideal(l) ,G ; 

j= psi (check) ; 
if (ecart (h)<ecart(j)) 
{ 

G[size(G)+l]=h ; 

} 

hl=h ; 
h=h-j ; 
kill s ; 

} 

} 

return (h) ; 

} 

SINGULAR Example 4.1 

ring r=0, (x,y) , Ds ; 

ideal i=x2,x4+x5+x6,x7,y2,y3+x8 ; 

poly f =x4y3+y5 ; 

WSNF(f, i) ; 

=> X5y3-x6y3-x8y2-xl2-xl3-xl4 
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ring r=0 , (x) , Is ; // example 3.6 
ideal i=x3+x4 ; 
poly g=x3+x6, x-x2; 
WSNFCg, i) ; 
=> 

2) Procedure to compute S-polynomials 

"sasbiSpoly procedure" : This procedure computes the generators of AR{G) (d 
fined in definition 3.10) which are S-polynomials. 

SINGULAR Procedure: 

LIB"elim. lib" ; // we need this library for "nselect" procedure 
proc sasbiSpolyCideal id) 

{ 

def bsr= basering ; 

ideal vars = maxideal(l) ; 

int n=nvars(bsr) ; 

int m=ncols(id) ; 

int z ; 

ideal p ; 

if (id==0) 

{ 

return (p) ; 

} 

else 

{ 

executeC'ring Rl=("+charstr(bsr)+") , (@y(l. .m) , "+varstr (bsr)+" 

(ds(in),ds(n));"); 
ideal id =imap(bsr , id) ; 
ideal A ; 

for (z=l; z<=in; z++) 
{ 

A[z]=lead(id[z])-(§y(z) ; 

} 

A=std(A) ; 

ideal kern=nselect(A,in+l,in+n) ; 
export kern, A ; 

setring bsr ; 
map phi= Rl , id ; 
p=simplify (phi (kern) , 1) ; 
return (p) ; 

} 



SINGULAR Example 4.2. 
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ring r=0, (x,y) , Ds ; 

ideal i=x2,x4+x5+x6,x7,y2,y3+x8 ; 

sasbiSpoly(i); 

[1] =x5+x6 

[2] =x8y 3+1/2x16 

3) SASBI BASIS construction algorithm 

" Sasbi procedure" : It is an iterative consequence of previous procedures to com- 
pute sasbi basis. 

SINGULAR Procedure: 

proc Sasbi (ideal id) 
{ 

ideal S.oldS.Red ; 
list L ; 
int z,n ; 
S=id ; 

while ( size(S) !=size(oldS)) 

{ 

L=sasbiSpoly(S) ; 

n=size(L) ; 

for (z=l; z<=n; z++) 

{ 

Red=L[l] [z] ; 
Red=WSNF(Red[l] ,S) ; 
oldS=S ; 
S=S+Red ; 

} 

} 

return (S) ; 

} 

SINGULAR Example 4.3. 

ring r=0, (x,y) , Ds ; 

ideal i=x2 , x4+x5+x6 , x7 , y2 , y3+x8 ; 

Sasbi (i) ; 

[l]=x2 

[2] =x4+x5+x6 
[3] =x7 
[4]=y2 
[5] =y3+x8 
[6] =x5+x6 
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